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ABSTRACT 

We explore possible systematic errors in the mass measurements of stellar mass black holes. We find 
that significant errors can arise from the assumption of zero or constant emission from the accretion 
flow, which is commonly used when determining orbital inclination by modelling ellipsoidal variations. 
For A0620-00, the system with the best available data, we show that typical data sets and analysis 
procedures can lead to systematic underestimates of the inclination by ten degrees or more. A careful 
examination of the available data for the 15 other X-ray transients with low-mass donors suggests 
that this effect may significantly reduce the black hole mass estimates in several other cases, most 
notably that of GRO J0422+32. With these revisions, our analysis of the black hole mass distribution 
in soft X-ray transients does not suggest any "mass gap" between the low end of the distribution and 
the maximum theoretical neutron star mass, as has been identified in previous studies. Nevertheless, 
we find that the mass distribution retains other previously identified characteristics, namely a peak 
around 8M Q , a paucity of sources with masses below 5M Q , and a sharp drop-off above 10M Q . 

Subject headings: black hole physics — X-rays: binaries 



1. INTRODUCTION 

Soft X-ray transients (SXTs) provide some of the 
strongest evidence for the existence of stellar mass black 
holes. In these systems, Eddington limited X-ray out- 
bursts occur over timescales of weeks to mont hs, followed 
by y ears to decades of X-ray quiescence (jChen et al.1 
1997). During quiescence, the optical light from the sys- 
tem is dominated by the companion star. Classic binary 
star analysis techniques can be used to determine the 
orbit of the companion star in quiescence, and thus the 
orbital parameters of the entire system, including prc- 
cise determin ations of the mass of the compact object 
(jOroszl [2003( 1. Many of these mass determinations are 
greater than the theoretical upper bound of neutron stars 
3Mq), and thus the compact objects are understood 
to be black holes. These black hole SXTs (BHSXTs) 
comprise most of the dynamically confirmed stellar mass 
black holes. 

Mass measurements of these systems have been used 
to explore the distribution of black hole masses in X-ray 
binar ies (|Bailvn et alj|1998t lOzel et al.|[2010t iFarr et all 
l2011h . In these studies, published values for the orbital 
parameters have been used as inputs for a Baycsian anal- 
ysis of the mass distribution of the compact objects. In 
all of these analyses, a significant mass gap between the 
maximum neutron star mass (3M©) and the low end of 
the black hole mass distribution (< 5M ) has been iden- 
tified. This is a curious result, since one would ordinarily 
expect that the mass distribution of black holes would 
be weighted toward the low end, as is the mass distri- 
bution of prc-supcrnova massive stars. Several theories 
have been proposed to explain the evolution of massive 
stars in binary systems and how the resulta nt supernova 
explosions might result in su ch a mass gap ()Brown et al.1 
120011 : iBelczvnski et al.ll20rl . 



Recently, however, the accuracy of the mass measure- 
ments for individual objects has implicitly been ques- 
tioned. Cantrell et al. (2010, hereafter C10) have re- 
analyzed all extant data for the prototypical BHSXT 
A0620-00. Prior to C10, published values for the or- 
bital inclination ra nged from i = 37°, im plying a black 



hole mass of 16M Q (jShahbaz et al.lll994h to i = 74°. im- 



plyin g a black hole mass of 4M Q (IFroning fc Robinson 



2001 ) with several inte rmediate values IjHaswell et al 



19931 : iGelino et alJ feOQl). C10 find that this wide range 



of incompatible results can be reconciled by more careful 
modelling of the ellipsoidal light curves of the companion 
star. In particular, C1 cull the data to inc lude only data 
in the "passive" state ([Cantrell et al.ll2C)08| ) , in which the 
lightcurves do not exhibit short-term non-ellipsoidal vari- 
ability. They fit the lightcurves with a model that in- 
cludes variable disk light and a hotspot that is allowed 
to vary in position and brightness. In this way C10 find a 
consistent value of the orbital inclination, i = 51.0±0.9°, 
which implies tobh = 6.6 ± 0.25M Q . 

It is becoming increasingly clear that other BHSXTs 
have active/passive state changes si milar to A0620-00, 
and m ay have variable hotspots (e.g.. iMa cDonal d et al.l 
1201 lh . However, the orbital parameter estimates in the 
literature generally do not take these effects into account. 
Thus it is possible that, like A0620-00, many of the BH- 
SXTs may have mass estimates that are inaccurate by 
considerable amounts. In this paper we examine the sys- 
tematic errors introduced by the presence of nonstcllar 
flux in the ellipsoidal lightcurves of BHSXTs. We find 
that there is likely a bias toward mass estimates that 
are higher than the true mass of the compact object. 
We find that the "mass gap" identified in previous work 
(jBailvn et al.iri998L [C^efet^l2?Tinl : IFarr et al.ll20lih on 
the mass distribution of BHSXTs can be accounted for 
by this systematic effect. 
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The outline of the paper is as follows. In §2 we sum- 
marize the primary sources of systematic uncertainty in 
black hole mass measurements. We quantify these un- 
certainties in §3 for A0620-00. In §4 we generalize the 
error estimates obtained for A0620-00 to other systems. 
§5 contains a re-evaluation of the mass estimates for 16 
BHSXTs. We use these revised mass estimates to ana- 
lyze the mass distribution in §6, and in §7 we conclude. 

2. SOURCES OF SYSTEMATIC ERROR IN MASS 
DETERMINATION 

The mass ttibh of the black hole in a BHSXT is de- 
termined by three parameters: the mass ratio of the sec- 
ondary star to the black hole, q = m*/mBH, the mass 
function /, and the orbital inclination i. Written in terms 
of these parameters, the black hole mass is 



sin i 

The error on the mass measurement thus depends on the 
error in /, q, and i. 

The mass function / generally contributes little to the 
systematic error on black hole mass. It is a function of 
the orbital period P and semi-amplitude of the radial 
velocity curve K, f = PK 3 / (2nG). For most cases, P is 
measured to sub-percent precision using photometric ob- 
servations that span many periods. The radial velocity 
curve amplitude is typically known to within 10%, with 
precision mainly limited by the resolution of the spec- 
troscopy. Neither measurement is strongly affected by 
systematic bias. A narrow Gaussian is thus a good ap- 
proximation of the error on /, which introduces a small 
random error on msH- 

The mass ratio q is inferred from the rotational broad- 
ening of spectral lines. By inspection of equation (1), 
we see that tobh is relatively insensitive to the mass ra- 
tio for q -C 1. In 11 of the 16 systems in our sample, 
q < 0.15 (see Table 2). For three systems with larger q, 
the mass ratio is well-constrained. The remaining two 
systems have only an upper limit, q < 0.5. Thus in all 
but these two systems the uncertainty in the measured 
mass ratio has a small impact on tobh- 

By far the largest source of systematic error is the or- 
bital inclination. Inclination is typically measured by an- 
alyzing ellipsoidal variability in the observed photometric 
lightcurve. The origin of this variability is gravitational 
distortion of the companion star. The star fills its Roche 
lobe, so as it orbits the black hole, the projected surface 
area and average temperature along the observer's line 
of sight is not constant. This gives rise to characteris- 
tic double peaks in the lightcurve, known as ellipsoidal 
variations, whose amplitude depends on the inclination. 
The largest amplitude occurs for systems edge-on to our 
line of sight (i = 90°), because that geometry maximizes 
the changes in projected surface area of the star with 
orbital phase. By contrast, no ellipsoidal variations can 
be detected for face-on systems (i = 0°). One can there- 
fore determine i by modeling ellipsoidal variability in the 
observed photometric lightcurve. 

The simplest model of ellipsoidal variability is just a 
simulated Roche-lobe filling star. This model, which we 
will call the 'star-only' model, is commonly employed 
in the literature under the assump tion that nonstellar 
sources of light arc negligible (e.g. iMartin et al.l 119951 : 



iGreene et all 120011: iGelino fc Harrison! l200l . It is not 
unusual, however, for nonstellar sources to contribute 
more than half the total flux from the system (e.g. 
IZurita et al.ll2002t lOrosz et alj|2004 C10). The ratio of 
nonstellar light to the total flux (hereafter denoted NSL 
fraction) is critical to measuring inclination accurately. 
Several sources contribute to the NSL fraction, including 
the accretion disk, hotspots on the disk, and potentially 
a jet. They can distort the shape of the photometric 
lightcurve in the following ways: 

1 . The accretion disk contributes a baseline flux which 
dilutes the amplitude of the ellipsoidal variations. 
To first order, increasing the disk flux is degenerate 
with lowering the inclination. 

2. The disk and jet can exhibit aperiodic changes in 
brightness (flickering) due to the accretion flow 
(C10). The timescale of this variability is short 
compared to the orbital period, so flickering can 
be misinterpreted as photometric error, particu- 
larly in folded lightcurves. However, flickering is 
not a white noise process, nor is it stationary, so 
binning and averaging the data will not reproduce 
the underlying ellipsoidal shape. 

3. The disk hotspot can also distort the lightcurve 
shape. The hotspot causes a peak in flux once per 
orbit, which leads to assymmctric lightcurves when 
superposed with double-peaked ellipsoidal variabil- 
ity. Morever, the position of the hotspot is not con- 
stant, so it can alternately increase or decrease the 
amplitude of the ellipsoidal variations for the same 
source (C10). 

Measuring inclination accurately is thus a delicate proce- 
dure that requires fitting a lightcurve with a consistent 
shape and a known NSL fraction, using a model that 
includes a disk and a hot spot. Using a less sophisti- 
cated model will result in systematically biased inclina- 
tion measurements. Moreover, because the lightcurve 
shape changes over time, different observing runs will 
yield different inclination measurements, even if the same 
model is used to fit the lightcurve. 

3. QUANTIFICATION OF SYSTEMATIC ERROR FOR 
A0620-00 

We now characterize the systematic error in inclination 
measurements due to the effects discussed in the preced- 
ing section. Our goal is to determine how inclination 
measurements obtained with a given model compare to 
the true inclination of a BHSXT. Specifically, we want 
to know the distribution of inclination measurements one 
would obtain from observing a source at many distinct 
times, and where the true inclination lies within that dis- 
tribution. This requires two pieces of information: one, 
an accurate inclination measurement, and two, a descrip- 
tion of the time variability of the source. 

Such an analysis is possible for the source A0620-00, 
one of the best-studied BHSXTs. Archival data spans the 
past 30 years, including eight observing runs and more 
than a decade of daily photometric monitoring from the 
Small and Moderate Apertu re Research Telescop e Sys- 
tem (SMARTS) consortium (jCantrell et al.l [20081 ). C10 
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use this extensive dataset to make an inclination mea- 
surement that accounts for the systematic errors dis- 
cussed in §2. They select 8 lightcurves that each maintain 
a consistent shape over the duration of the observations. 
They fit every lightcurve with a 11-parameter model that 
includes an accretion disk with variable temperature, size 
and flaring angle and a hotspot with variable tempera- 
ture, size, and position. During the fitting procedure, 
they constrain the disk to contribute the spectroscopi- 
cally determined NSL fraction. The best-fit inclinations 
to all 8 lightcurves are statistically self-consistent, giving 
a weighted average of i = 51.0 ± 0.9°. We assume this 
value is unbiased by systematic error. 

Given an accurate inclination measurement for A0620- 

00, we now determine how inclination estimates vary 
with time. First we select all subsets of existing data that 
resemble plausible observing runs. We fit each of these 
lightcurves with a star-only model to obtain a sample of 
inclination estimates. We chose a star-only model both 
because it is commonly used and because it introduces a 
straightforward systematic error. 

3.1. Distribution of Inclination Measurements from 
Archival A0620-00 Data 

For our analysis, we use data collected by CIO. We 
select all subsets of this data that could be obtained on a 
typical observing run, i.e., many observations in a short 
period of time. Specifically, we require each subset to 
contain at least 20 data points over 7 days, with gaps in 
phase no larger than 0.1. We also specify that subsets 
do not overlap, ensuring that each point is only counted 
once. We restrict the sample to lightcurves of compara- 
ble quality to previously published BHSXT lightcurves 
by binning the data in 30 phase bins and removing all 
lightcurves with average bin deviation greater than 0.03 
mag. Our final sample consists of 57 lightcurves taken 
with four different filters: V, I, H, and W (a wide band- 
pass centered at 4700 A). 

We fit each lightcurve with a star-only model to ob- 
tain a distribution of inclination estimates. To fit the 
data, we use the E c lipsin g Lightcurve Code (ELC) of 
lOrosz fc Hauschildtl (|2000D . Our model star has a mean 
temperature of 4600 K and a gravity darkening exponent 
of 0.10, consis tent with a K5 spectra l type and a convec- 
tive envelope (jLucvll 19671 : lGravlfT9 92T) . Limb darkening is 
computed directly from the model atmosphere. For each 
lightcurve in our dataset, we obtain an inclination esti- 
mate i from the model that gives the lowest Xred- The 
distribution of i for all the lightcurves is shown in Figure 

1. Note that we distinguish the estimated inclination i 
from the true inclination, denoted i. We find i ranges 
from 33.9° to 52.4°. Over 90% of the measurements fall 
below the C10 measurement of i = 51°. Such systematic 
underestimation is expected: nonstellar flux almost al- 
ways dilutes the amplitude of the ellipsoidal variations, 
mimicking a lower inclination. The few lightcurves with 
i > 51° may be the result of constructive interference 
between the hotspot and the ellipsoidal variability. The 
distribution of i is bimodal, with one peak near 40° and 
the other near 50°. We discuss potential mechanisms for 
this bimodality in the next section. 

3.2. Distribution of Inclinations from Simulated 
A0620-00 Lightcurves 
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Fig. 1. — Histogram of inclination estimates obtained by fitting 
a star-only model to extant A0620-00 lightcurves. Note that 90% 
of the probability mass lies below the C10 inclination measurement 
i = 51°, marked in this plot with a thick dashed line. 

The bimodality in the distribution of inclinations seen 
in Figure 1 is an intriguing result because it illustrates the 
potential impact of nonstellar flux on inclination mea- 
surements. To explore how different factors contribute 
to the shape of the distribution, we introduce a method 
to simulate A0620-00 lightc urves. 

We build on the work of ICantrell et al.l (|2008ft . who 
analyze the time variability of A0620-00 using long- 
term photometric monitoring from the SMARTS consor- 
tium. They identify two distinct states — labeled "ac- 
tive" and "passive" — in the quiescent, optical A6020-00 
lightcurve. In the passive state there is minimal aperiodic 
variability, so the lightcurve shape is stable from night to 
night. By contrast, flickering is much more pronounced 
in the active state: active lightcurves are brighter, bluer, 
and more variable, possibly due to increased accretion 
activity. States persist for several months, but transi- 
tions from state to state occur on sub-night timcscales. 
A0620-00 is active roughly 70% of the time (C10). Be- 
cause active and passive data behave differently, we sim- 
ulate the two states separately. 

3.2.1. Simulation of Passive Data 

Although passive-state lightcurves do not change shape 
on short timescalcs, passive states separated by a period 
of activity do not necessarily have the same lightcurve 
shape. The 8 lightcurves fit by C10 were passive, but 
they had significantly different shapes from each other. 
The main difference between the model fits were the 
hotspot parameters, suggesting that shape changes in 
passive lightcurves are due to the changing temperature 
and position of the hotspot. Therefore, to model passive 
data, we select hotspot parameters from the range of fits 
in C10, listed in Table 1 (J. Orosz, private communica- 
tion). We find no correlation between any of the parame- 
ters, so we choose each parameter uniformly between the 
minimum and maximum value given in Table 1. We use 
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Fig. 2. — Simulated active V-band lightcurves (points). The 
best fitting star-only models (lines) have inclinations i = 30.4°, 
37.4°, and 44.9°, with larger ellipsoidal variability corresponding 
to higher inclination. Each lightcurve consists of a single night of 
simulated data. 

ELC to generate 500 passive V- and H- band lightcurves 
with the specified hotspot parameters. We then add pho- 
tometric errors drawn from a normal distribution with 
mean and standard deviation 0.03 mag. A representa- 
tive sample of simulated lightcurves is shown in Figure 2. 
The lightcurves are fit with a star-only model with pa- 
rameters described in §3.1. All fits had Xred < 3.0. We 
consider the quality of these fits acceptable, as Xr C d ^ as 
commonly ex ceeded 2.0 in previously published inclina- 
tion fits (e.g. IShahbaz et alj|1996t iBeekman et al.|[l997l : 
Ivan der Hooft et al.l I1998D The distribution of best-fit 
inclinations for each filter is shown in the top panel of 
Figure 3. These distributions reveal that the A0620-00 
hotspot has a greater effect on the lightcurve shape at 
longer wavelengths. For both filters, the distribution is 
peaked around i = 51° with some left skew. For an esti- 
mator i obtained by fitting a star-only model to passive 
data, the true inclination i is given by i = rtj '4 for H- 
band data and i = for V-band. The limits denote 
68% confidence. These distributions reproduce the sec- 
ond mode in Figure 1, confirming the suggestion that the 
hotspot can increase or decrease the amplitude of the el- 
lipsoidal variations. The distributions may be artificially 
narrow: the range of hotspot parameters listed in Ta- 
ble 1 is derived from only 8 fits, which probably do not 
span the entire parameter space. In addition to possible 
underestimation of the hotspot variability, we neglect all 
other disk and jet variation. We therefore expect the 
distribution of 1 for passive lightcurves to be somewhat 
broader than the results quoted above. 

3.2.2. Simulation of Active Data 

For active-state data, we model the stellar and non- 
stellar components of the flux separately. We simulate 



Fig. 3. — Kernel density estimates for inclination obtained from 
fitting a star-only model to passive and active lightcurves (top and 
bottom panels, respectively). Solid lines correspond to simulated 
H-band data and dashed lines to simulated V-band data. All four 
samples contain 500 lightcurves. The vertical gray line marks the 
inclination measurement obtained by C10. 



TABLE 1 
Hotspot parameters 



Fit 


^spot 


^spot 


r cu t 


Wspot 


1 


7.14 


341.7 


0.37 


33.9 


2 


2.53 


200.8 


0.27 


31.7 


3 


1.53 


19.3 


0.87 


90.0 


4 


2.41 


199.4 


0.89 


89.7 


5 


5.70 


190.3 


0.07 


70.7 


6 


4.20 


102.5 


0.66 


85.9 


7 


6.36 


72.5 


0.54 


39.4 


8 


5.50 


37.7 


0.38 


20.5 



Notes: This table lists the hotspot parameters for the ELC fits 
in C10. The ratio of hotspot temperature to disk temperature 
is Sspot- The hotspot is centered at azimuthal angle 6 Bpo t and 
extends from the outer disk to an inner radius r cu t . The angular 
radius is given by ui spo t. 

the stellar component with ELC using the same temper- 
ature and gravity darkening parameters as in §3.1. We 
set the average stellar magnitude equal to the zero-disk 
magnitude found in C10. 

Next we add a nonstcllar component to the stellar 
lightcurve. We model the nonstellar flux with a bro- 
ken power law power spectral density (PSD). The PSD 
cannot be determined directly from existing lightcurves 
due to severe daily and yearly aliasing. We therefore 
develop an alternative method to characterize the time 
variability of the nonstellar flux. We then determine the 
broken power law parameters that best reproduce that 
variability. 
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We begin by characterizing the time variability of the 
nonstellar flux. Our data consist of H- and V-band 
active-state lightcurves from the Small and Moderate 
Aperture Research Telescope System (SMARTS) c onsor- 
tium, originally published in ICantrell et al.1 (j2008t ). The 
H-band dataset contains 1109 points over 2227 days and 
the V-band dataset contains 634 points over 2755 days. 
We subtract the ellipsoidal variability from this data 
using the stellar lightcurves produced with ELC. The 
median time between observations is one day; however, 
there are large seasonal gaps in the data when the object 
is obscured by the sun. Even when it is visible, sam- 
pling is somewhat sporadic. Due to such uneven sam- 
pling, conventional methods to characterize time vari- 
ability fail. We therefore use a modification of the binned 
autocorrelation function, denoted Ca(i"), to describe the 
time variability. 

For the time series {(ii, ft), (t n , /„)}: 



n 



(2) 



where 



• t is the lag, 

• /(o ± b) is the average of the set {fj} such that 
{tj} is within the range a ± b 

• and N is the number of sets where {fj} ^ 0. 

In other words, Ca(t) is the average squared difference 
for points separated by r ± A. If there are no points 
separated by r ± A (i.e., N = 0), Ca(t) is undefined. 
To reduce the number of pairs that are likely to be cor- 
related, we do not include any (U, fi) more than once 
for a given lag r . For example, if r = 10.0 ± 0.05 and 
t = {0.0, 0.01, 10.0}, we remove the pair {0.01, 10.0} be- 
cause t = 10.0 has already been included. The error on 
Ca for a given r is the sample standard deviation for the 
N squared differences. Note that these errors are not 
Gaussian, because pairs of points separated by r ± A are 
not statistically independent. 

The SMARTS observations provide sufficient coverage 
to constrain Ca(t) for lags on the order of 10 _1 to 10 3 
days. The median time difference between observations 
for both H- and V-band data is 1.0 ± 0.05 days, so we 
evaluate Ca(t) starting at r = 1.0 and increasing by 
factors of 2.0 to r = 1024.0. We choose A = 0.1 days. 
Variability on timescales < 0.1 days is small compared 
to variability for r > 1.0, so this choice of A does not 
bias the calculation of Ca- For r < 128.0 days, this A 
gives N > 300 for both filters. On longer timescales, 
there are fewer pairs of points from which to choose, so 
we exclude r with N < 15. To constrain variability on 
short timescales, we also compute Ca(t = 0.1 days) with 
A = 0.01, obtaining N — 25 for the V-band data and 
N = 31 for H. The results of the calculation of Ca are 
shown in Figure 4. 

We next determine the broken power law PSD that 
best matches the Ca statistic for each data set. We as- 
sume a PSD S(f) of the form: 



S(f) 



1 

7° 



if / < h 
if / > fb- 



(3) 



where / is frequency, fb is the break frequency, a is 
the slope at frequencies less than fb, and /3 is slope at 
frequencies greater than ft,. 

To find the power law parameters that best repro- 
duce the observed nonstellar lightcurve, we implement 
a grid-based x 2 minimization routine. We loop over a, 
j3, and fb and simulate a lightcurv e from each set of 
param eters following the method of iTimmer fc Koenid 
(|1995h . We normalize the lightcurve so that it has the 
same 25% and 75% quartiles as the observed nonstel- 
lar lightcurve. This normalization ensures that outly- 
ing data points do not bias the amplitude of the sim- 
ulated lightcurve. We then compute Ca(t) for the 
simulated data and calculate its goodness of fit to Ca 
for the observed data. We find the power law pa- 
rameters that give the mininum x 2 are ( a i P\ fb) = 
(-0.8tg;J, -1.5±g;£, 1.2±g;£ x 10~ 8 Hz) for V-band data, 
with x ?ed = 1-7, and'(-0.8t°;£,-1.6±°j,5.0±j*;° x 
10~ 9 Hz) for H-band, with x 2 cd = 2.7. The results for the 
observed and best-fit simulated lightcurves are shown in 
Figure 4. We obtain the uncertainties by holding two 
parameters fixed and determining the range of the third 
parameter over which x 2 increases by a factor of two. 
These ranges do not represent confidence intervals, nor 
do we expect x? c d ~ 1 f° r cither data set, because the 
errors on Ca(t) are not Gaussian. Nevertheless, varying 
the parameters within the quoted ranges does not signif- 
icantly change the distribution of i obtained from fitting 
the inclination of the simulated lightcurves. 

To obtain mock active-state data, we simulate fifteen 
2500-day long lightcurves in V and H using the best-fit 
power law parameters. The lightcurves are sampled at 
10-minute intervals. We normalize each lightcurve such 
that the 25% and 75% flux quartiles match the observed 
SMARTS data. We limit the length of the lightcurves 
to 2500 days because Ca(t) is poorly constrained for 
t > 1000. In addition, the normalization should be de- 
termined using lightcurves of comparable length. After 
normalizing the nonstellar lightcurves, we add back the 
stellar component of the flux. 

In order to select sample data from the 2500-day long 
simulated lightcurves, we recreate plausible observing 
runs in the following way: first, we remove data points 
during "daylight" hours, assuming 8 hours of viable 
observing time per 24-hour period. We simulate bad 
weather by removing 4-hours segments of data with prob- 
ability 20%. Samples of the lightcurves are then selected 
with with mock observing runs separated by 30-day in- 
tervals. Each run has equal probability of being one, two, 
or three consecutive days long. We exclude the run from 
the final sample if it has gaps in phase greater than 0.1. 
We also bin each run in 30 phase bins and restrict the 
sample to those runs with average bin standard devia- 
tion less than 0.03 mag and none greater than 0.1 mag, 
as we did for subsets of real data. We select 500 of the 
remaining lightcurves for each filter such that the sam- 
ples are the same size. Several representative lightcurves 
are shown in Figure 2. We fit each of the lightcurves in 
the sample with the star-only model described in §3.1. 
The resulting estimates of inclination, again denoted i, 
are shown in the bottom panel of Figure 3. The distri- 
butions are approximately normal. We fit a Gaussian 
to each distribution and find i ~ 7V(38.2°, 4.0°) for the 
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Fig. 4.— The statistic Ca for V- and H-band SMARTS data 
(points in the top and bottom panel, respectively) and data gen- 
erated from the best-fit power law (lines). The value of Ca repre- 
sents the average squared difference in flux between data points 
separated by time r ± A. The error bars represent the stan- 
dard deviation of Ca for the N points associated with each r. 
The choice of r for lags > 256.0 is not identical for the two 
data sets because they are sampled slightly differently and we re- 
quire N > 15 for each r. The best-fit power law has parame- 
ters (a, 0, f b ) = (-0.8,-1.5,1.2 X 10~ 8 Hz) for V-band data and 
(-0.8,-1.6,5.0 x 10" 9 Hz) for H-band. 

V-band fits and i ~ 7V(43.2°, 3.9°) for H-band, where 
the notation iV(/z, a) indicates a normal distribution with 
mean [i and standard deviation a . The relation between 
true inclination i and the estimator i is then given by 
i ~ N(i + 12.8°, 4.0°) for V and i ~ N{% + 7.8°, 3.9°) for 
H. 

3.3. Lessons from A0620-00 

Using simulated data, we find that fitting a star-only 
model to passive state lightcurves results in unbiased in- 
clination estimates, approximately normally distributed 
around the true inclination. The scatter can be explained 
by a hotspot with changing position and temperature. 
On the other hand, fitting active data introduces a bias 
towards artificially low inclination estimates. The typi- 
cal underestimation is around 8° for IR lightcurves and 
13° for optical lightcurves. The source of the bias is a 
significant NSL fraction that increases with shorter wave- 
length. 

These results reproduce the bimodality in inclination 
estimates obtained from past observations (shown in Fig- 
ure 1). The optical state of the source is unknown 
for many of these observations, but the first mode near 
i ~ 40° is comparable to the fits to active data and the 
second mode around i ~ 50° is consistent with the pas- 
sive fits. In addition, the relative weight of the two modes 
(roughly 2/3 of the probability mass centered on the first 
mode) matches the C10 observation that A0620-00 is ac- 
tive around 70% of the time. 



4. GENERALIZATION OF SYSTEMATIC EFFECTS TO 
OTHER SYSTEMS 

We will now use the description of systematic effects 
for A0620-00 as a framework to estimate the systematic 
error on inclination for other sources. We treat passive 
and active data separately. 

For passive A0620-00 data, we found that inclination 
estimates i obtained with a star-only model are related 
to the true inclination i by i = i^ 1 ■ We assume a simi- 
lar relation is valid for other systems, because true stellar 
ellipsoidal variability is dominant in passive lightcurves. 
Thus we expect inclination measurements made using 
passive data to be unbiased (centered on the true in- 
clination) for all systems. We choose o\ = o-i = 3.0°. 
The a values we obtained for A0620-00 are not symmet- 
ric and are somewhat smaller; however, as discussed in 
§3.2.1, they are most likely underestimates, so a = 3.0° 
is a conservative approximation of the error. 

For active A0620-00 data, we found 

i~7V(? + £,<r). (4) 

To generalize this result to other sy terns, we scale the 
values of £ and a based on the system's orbital param- 
eters. We find that £ can be computed from the NSL 
fraction, which we denote <f>. We discuss this approxima- 
tion in §4.1 and describe a method to estimate the NSL 
fraction in §4.2. We chose a to scale linearly with £. 

4.1. Dependence of Inclination Measurements on the 
NSL Fraction 

We expect that the bias in inclination measurements 
for active data is determined primarily by the NSL frac- 
tion, and that flickering and hotspots are secondary 
effects. To demonstrate this, we simulate a set of 
lightcurves with two components: a star and a constant 
offset flux representing the NSL fraction. We choose 
i = 51° for the stellar lightcurve and NSL fractions rang- 
ing from 0.0 to 0.9. We fit these diluted lightcurves with 
a star-only model to obtain i as a function of the NSL 
fraction, <f>. A third-order polynomial fit to the results is 

i- 1 = 0.2 + 28.60 - 15. 30 2 + 27.O0 3 . (5) 

where we take i = 51°. We plot this fit in Figure 5 and 
overplot the points (</>, i) for each of the simulated A0620- 
00 active lightcurves. These points fit the curve closely: 
the simulated active data have i within ±2.7° of the value 
obtained from adding a constant NSL fraction (limits 
denote 68% confidence). We therefore conclude that the 
NSL fraction is indeed the dominant factor determining 
the bias £ = i — i, with some scatter introduced by the 
disk hotspot and flickering. 

4.2. Estimating the NSL Fraction Using Orbital 
Parameters 

Given that the NSL fraction determines the bias £ to 
first order, we would like to estimate the typical NSL 
fraction for sources other than A0620-00. Qualitatively, 
we expect a lower NSL fraction for systems with rela- 
tively hotter stars and a higher NSL fraction for systems 
with relatively larger accretion disks. To quantify these 
relationships, we focus on three observable parameters: 

1. Spectral type of the secondary star. The hotter the 
star, the greater its total flux contribution and the 
lower the NSL fraction. 
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Fig. 5. — Inclination estimates from a star-only model as a 
function of phasc-avcragcd NSL fraction for simulated active H- 
and V-band lightcurves (open and closed circles) . The line gives the 
inclinations fit to stellar lightcurves with a constant NSL fraction 
added. The active lightcurves have i within ±2.7° of the black 
curve at 68% confidence. The dashed line marks the inclination 
measurement of C10, i = 51.0°. 



2. Mass ratio. The mass ratio determines the relative 
size of the Roche lobes in the system. Assuming 
the fraction of the Roche lobe occupied by the disk 
is constant over all systems, a system with higher 
q has a relatively larger disk and correspondingly 
higher <fr. 

3. Inclination. As i increases, the projected area of 
the disk decreases, thus lowering <j). 

To scale the NSL fraction based on the system param- 
eters listed above, we use the Stefan-Boltzmann law. For 
an optically thick source, the flux scales with the area of 
the emitting region and the temperature to the fourth 
power. For convenience, we work with the ratio of stellar 
to nonstellar flux, which we denote p. This parameter is 
related to the NSL fraction by p = l/cf) — 1. We expect 
p to scale as: 



poc 



■ rp \ 4 / p \ 2 
star \ / -^star 




sec i 



sec i, 



(6) 
(7) 



where T and R denote temperature and effective radius, 
with subscripts denoting the star and disk. Equation (7) 
repla ces the ratio of radi i with q using an approximation 
from iFrank et al.l (|2002t ) . The sec i factor accounts for 
the orientation of the accretion disk, assuming it lies en- 
tirely in the orbital plane. We approximate that the star 
is spherical, so its flux is independent of inclination. 

We can use the expression in equation (7) to estimate 
the NSL fraction of any source, provided we know the 



constant of proportionality. To solve for the constant, we 
use the known parameters for A0620-00. We use the esti- 
mates of q, i, and T sta r from C10. We estimate p from the 
simulated active A0620-00 data, obtaining p = 1.0 ±0.25 
in the optical and 2.3 ± 0.6 in the IR. The remaining 
unknown is Tdisk, which we choose to absorb into the 
constant of proportionality. By neglecting the Td; s k de- 
pendence, we effectively assume that all sources have the 
same disk temperature as A0620-00. This approxima- 
tion is reasonable because the internal dynamics of the 
disk determine its temperature, not the properties of the 
gas as it passes through the inner Lagrange point. In 
addition, this assumption is substantiated by empirical 
evidence from three cases: the nonste llar light is redder 
than the B star in SAX J1819.3-2525 (iMacDonald et al.l 
120111 ). bluer than the K star in A0620-00 (C10), and 
approximately constant across the optical spectrum for 
the F star in GRO J1655-40 (Cantrell, private commu- 
nication). Under this assumption, we can solve for the 
proportionality constant in (7) and thus infer p seci for 
any system with known q and T s t a r- 

We can use our estimate of p sec i to infer the bias £ in 
a star-only inclination estimate, i. First, we generate a 
set of stellar lightcurves and add a constant flux offset to 
each such that the ratio of stellar to nonstellar flux equals 
p seci. We vary i from 0° to 90° in 1° intervals. Second, 
we fit each lightcurve with a star-only model. We then 
find the lightcurve whose inclination estimate is closest 
to i. Because we know i for this lightcurve, we can obtain 
the bias from £ = i — l and the NSL fraction according to 
equation (5). To evaluate the accuracy of this procedure, 
we compute the NSL fraction from pseci for all sources 
with star-only inclination estimates. We compare these 
values to spectroscopically determined NSL fractions in 
Figure (6). There is good agreement between our predic- 
tions and the observed values, with 7 of 11 predictions 
within lcr from the observed value. 

The spread in the NSL fraction of A0620-00 introduces 
a small amount of scatter in our estimate of £, generally 
< 2°. This error is unrcalistically narrow because it does 
not include any uncertainty due to other components of 
the lightcurve, such as flickering. To estimate the a in 
equation (4), we choose the linear relation a = 3.0 + 
0.115£. This scaling reduces to the expression for passive 
data in the case £ = 0.0° and ensures that a = 3.9° for 
£ = 7.8° (i.e., the values we find for simulated active 
H-band lightcurves). 

4.3. Additional Constraints on Inclination 

There are two additional constraints on inclination we 
can obtain. First, we set an upper limit based on the 
absence of eclipses. The visibility of eclipses depends on 
the size of the accretion disk and the mass ratio, q. We 
assume a conservative disk radius of 0.5i?2 and use ELC 
to determine the eclipse limit as a function of q. Systems 
with lower q eclipse at higher inclinations, so we use low- 
end estimates of q to establish secure upper limits on i for 
the all BHSXTs with measured mass ratios. Specifically, 
we use the minimum value for mass ratios described by 
a uniform distribution and the value two standard de- 
viations below the mean for mass ratios with a normal 
distribution. The upper limits calculated according to 
this procedure are listed in Table 2. 

The second constraint on inclination is a lower limit ob- 
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Fig. 6. — Observed versus predicted optical NSL fractions 
(points). The gray line is a visual aid resprescnting perfect agree- 
ment between the observations and predict ions. The numbers indi- 
cate the reference for th e observation: ( 1) [O rosz ct al. ( 19 93), (2) 
" ||2000|'), (3) ICasares eTaLl lft995l), (4) IShahbaz ct al. 
Casares et al.l 111993). (6) IFilipp ciiko ct al. ( 199 9|), (7) 



IWebb et al 



(199 9]), (5) 



IMacDonald et al.l i20TW, (8) IQrosz et al.l (120021), (9 ) IQrosz et al.l 
(120041) , (10) IQrosz et all 11999) . (11) IWagner etall 1 120011) . The 

vertical error bars indicate either la error or an allowed range. 
The horizontal error bars reflect the uncertainty due to the vari- 
ation in the NSL fraction of A0620-00, on which the predictions 
are based. The excess in the observed NSL fraction for (9) may 
because the source (GRO J1650-500) was not in full quiescence dur- 
ing the observation. The lack of agreement for source (11) (XTE 
J1118+480) may be due to this system's high inclination; for large 
i, our assumption that the accretion disk lies entirely in the orbital 
plane may be invalid. 



taincd by assuming the secondary star has a mass equal 
to or less than that of a main sequence star of its spectral 
type. A Roche-lobe filling star is out of thermal equi- 
lbrium, so the relationship between mass, radius, and 
surface temperature may be quite different from that of 
spherical stars. Thus far, undcrmassive secondaries have 
bee n observed in several sour ces, including GRO J1655- 
40 (|van der Hooft et al.lfl997h and A0620-00 (CIO). We 
therefore interpret the standard mass for a secondary 
star's spectral type as an upper limit to its true mass. If 
q is known and m* is overestimated, tobh is also over- 
estimated. By inspection of equation (1), we find that 
overestimating tobh results in an underestimation of i for 
fixed q and /. There is some uncertainty in the spectral 
type for most of the sources in our sample, so to obtain 
the most stringent lower limit, we calculate i using the 
mass of the star with the earliest allowed spectral type, 
along with the low-end estimates for q and / described in 
the preceding paragraph. If the low-end estimate for q is 
less than 0.0, the lower limit on i is also 0.0°. The lower 
limits on inclination for all sources are listed in Table 2. 

5. RE-EVALUATION OF INCLINATION AND MASS 
ESTIMATES FOR TRANSIENT BLACK HOLE BINARIES 



Using the procedure established in the previous sec- 
tion, we examine the inclination measurements for the 
black holes in our sample. Much of the previous work 
on black hole X-ray binaries has had the goal of estab- 
lishing a lower limit on the mass of the compact object 
to demonstrate that it is a black hole and not a neutron 
star. As a result, the conservative lower limits on mass 
are often secure, but the masses and their errors may not 
be accurate. There a number of standard practices that 
may lead to erroneous measurements, which we briefly 
outline before discussing individual objects in detail. 

One of the most frequent assumptions when fitting 
inclination is that the non stellar flux is neglibl e , par- 
ticularly in the infrared (IShahbaz et al.l [l99l 119961: 
Beekman et~aTl Il997t IGelino et al.l 120011: iGreene et al l 



20011 : iGelino fc Harrison! I2003D . In many cases, this as- 



sumption has been justified by claiming the disk spec 
trum can be modeled as a po wer law wit h negative 
slope, based on the precedent of lOkd (|1977l ). The ac- 
cretion disk thus contributes less total flux in the IR 
than in the optical. However, this approach has been 
called into question by recent measure ments of signifi- 
cant IR NSL fractions i n some sources (| Reynolds et al.1 
l2008t IGelino et ai]|2010D . In addition, there is evidence 
that the NSL fraction increases with wavelength for some 
objects (jOrosz fc Bailvnlll997HRevnolds et al.ll2007l ). In- 
clinations obtained from fitting a star-only model should 
therefore be treated with caution. 

Another common practice when fitting inclination 
is to combine l ightcurves from two or more observ- 
ing runs (viz. IShahbaz et al.l 119941: iRemillard et all 



1996 a|: ["Gelino fc HarrisonT 20031: lloannou et al 



Orosz et a l. 2004; Gelino et alll2006l : ICasares et al.l 



2004; 



2009; 



Orosz et "all 1201 ll ). However, due to aperiodic flicker- 



ing and migratory hotspots, it is unlikely that combin- 
ing and/or binning data reproduces the true underly- 
ing lightcurve shape. In addition, if the datasets be- 
ing combined have non-overlapping phase coverage, large 
changes in the lightcurve shape may not be detectable. 
Binning data that has significant intrinsic aperiodic vari- 
ability tends to depress the amplitude of the ellipsoidal 
variations and thus lead to systematically low inclination 
measurements. 

A third potential source of error is using a NSL frac- 
tion determined at a different time from the photomet- 
ric observations. The nonstellar flux can exhibit signif- 
icant variability on timescales of months, as we show 
in §3.3 for A0620-00. Similar variability has been ob- 
served in many other sources, including XTE J1550-564, 
SAX J1819.3-2525 GRO J0422+32, GS 1354-64, and 
XTE J1859+226 ([Reynolds et all 120071: ICasares et all 



2009tlCorral-Santana et alj2011t IMacDonald et al.ll2011fc 



Orosz et al.ll2011| ). Given these changes in brightness, it 
is useful to determine the NSL fraction simultaneously 
with photometry to obtain an accurate inclination mea- 
surement. 

In the following section, we examine the data available 
for each of the 16 BHSXTs in our sample. In partic- 
ular, we attempt to determine whether the source was 
observed in an active or passive state. We then deter- 
mine an appropriate value and error on the inclination, 
as described in §4. In some cases, this results in values of 
i and hence of the black hole mass that are quite different 
from those in the literature. More extensive observations 
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will be required to check and refine these estimates. 

5.1. 4U 1543-47 

4U 1543-47 is one of the few source s with near simulta - 
neous photometry and spectroscopy. lOrosz et al.1 ()1998l ) 
analyze B-, V-, and I-band photometry from 1998 June 
28 - July 4 and spectroscopy from 1998 July 1-4. The 
authors determine the NSL fraction in B, V, R, and I 
from the observed spectrum and use the measurements 
as constraints when fitting models of ellipsoidal variabil- 
ity to the lightcurve. Due to the relative brightness of 
the A-star secondary, we expect the NSL fraction to be 
small. A larger source of error is the mass ratio, which 
the authors include as a free parameter in the model. 
They find 24° < i < 36° at the 3er level, but note the 
possibility of additional systematic effects. We therefore 
use this range as the boundary of a uniform distribu- 
tion. A more precise measurement, i ~ 7V (20.7°, 1.5°), 
was r eported in conference proceedings by lOrosz et al.l 
(HH). We do not use this measurement because there 
is no published record of the lightcurve. However, we 
note that using this smaller inclination would increase 
the most likely black hole mass by more than a factor of 
two. 

5.2. A0620-00 

We adopt an inclination measurement of 51 ± 0.9° for 
A0620-00, based on the analysis of CIO. See §3 for a 
discussion. 

5.3. GRO J1655-40 

GRO J1655-40 is a well-stu died source, wit h very 
regular quiescent lightcurves (lOrosz fc Bailvnl H997t 
iGreene et all 120011: iBeer fc Podsiadlowskil 1200a) . The 
consistency in lightcurve shape exemplifies the connec- 
tion between the stellar temperature and the NSL frac- 
tion that was discussed in §3.3: GRO J1655-50 has a 
bright F star secondary, so we expect that it has rel- 
atively small NSL fraction, and consequently, a stable 
lightcurve shape. Detailed studies of optical/IR emis- 
sion and orbital parameters of G RO J1655-40 in full qu i- 



69.0° with Xr 



1.6. Recent unpublished mea- 



escence have been p ublis hed by lOrosz fc Bailvn 
IGreene et all (l200lh a nd IBeer fc Podsiadlowski 



1997 ) 



2002) 



The lOrosz fc Bailvnl (jl997t ) study gives a tight con- 
straint on i because their modeling requires a partial 
eclipse of the secondary to account for the depth of the 
minim um at phase 0.5. By contrast, the IGreene et al.l 
(|2001[) analysis uses a more sophisticated model for limb- 
darkening that all ows for a deep primar y minimum with- 
out eclipses. The IGreene et al.1 ()2001l ) analysis use op- 
tical and IR lightcurves collected between 1999 July 7 
and October 30. Over the course of observations, the 
lightcurves remained stable within the range of photo- 
metric error, a few hundredths of a magnitude. The 
lightcurves are also cons istent with those repor ted from 
several years before by lOrosz fc Bailvnl (|1997|) . Thus 
these lightcurves appear to be passive. The data are 
fit with a star-only model , giving a best-fit incl i nation 
i = 70.2° with Xr ed = 1.6. IBeer fc Podsiadlowskil (120021) 
obtain a similar result by reanalyzing the lOrosz fc Bailvnl 
(119971 ) data. They use a model that includes an accretion 
disk and distance constraints obtained from the kine- 
matics of the radio jet, finding a best fit inclination of 



surements of the NSL fraction (Cantrell, private com- 
municati on) suggest that the distance to the source ob- 
tained bv IBeer fc Podsiadlowskil (|2002T) is more accurate 
tha n that of IGreene et all (|200lf) . We therefore adopt 
the IBeer fc Pod siadlowski (200 3) m e asurem ent , which 
is compatible with the IGreene et all (|2001[) result, and 
assume the source was passive during the observations, 
thus taking i ~ AT(69.0°, 3.0°). 

5.4. GRO J0422+32 

GRO J0422+32 has an unusually late-type secondary 
(jFilippenko et al.1ll995f ). and thus we expect the disk con- 
tribution to be especially strong. The inclination has 
be en frequently discussed since it was first measured 
by lOrosz fc Bailvnl ()1995l ). That work finds i > 45° 
based on I-band photometry from 1994 October 27-28. 
The lightcurve exhibits noticeable shape changes be- 
tween the two nights and a difference in mean I mag- 
nitude of 0.05. This type of variability is typical of the 
active state. A lower limit on i is obtained by fitting 
a star-only model to ellipsoidal variability with ampli- 
tude 0.15 mag, an approximate amplitude determined 
by considering each night in the observed lightcurve sep- 
ar ately. This measurement is consistent with the work 
of iGehno fc Harrison! (|2003t) . who find i = 45 ± 2° from 
fitting a star-only model to J, H and K' lightcurves. 
There is evi dence, however, for sub stantial IR disk con- 
tamination: iRevnolds et al.l (|2007t ) detect no ellipsoidal 
variability in H and K lightc urves that have mea n mag - 
nitudes consistent with the iGelino fc Harrison! ()2003f ) 
ligh tcurves. These con fl icting results su ggest that both 
the IGelino fc Harrison! (!2003[ ) and the IRevnolds et al.l 
(|2007f ) li ghtcurves are in th e activ e state, so we inter- 
pret the IGelino fc Harrison! (|2003l) measurement as a 
lo wer limit. These r esults are consistent with the work 
of iFilippenko et al.l (jl995D . who find i = 48 ± 3°, as- 
suming a normal mass M2 V secondary and mass ratio 
q = 0.1093 ± 0.0086. As discussed in §4.3, assuming a 
normal mass secondary gives a lower limit on inclination. 

In contrast to the above me asurements, there are sev- 
eral works which find i < 45°: ICasares et al.l (|1995l ) ob- 
tain % = 30 ± 6° for q = 0.1 j using I-band data and as- 
suming zero disk contribution. However, there are signif- 
icant gaps in phase coverage near the primary maximum 
at phase 0.25. Anoth er low measurement is obtained by 
iCallanan et al.l (|1996t) . who find i < 45° for an assumed 
disk contamination less than 0.2. This NSL fraction 
is based on what is considered typical for other SXTs, 
which may not be applicable to this source because of 
the late-type secondary. Moreover, the lightcurve they 
fit consists of 14 nights of data between 1994 September 
11 and 1995 January 9 that are combined and binned. 
Given the substantial aperiodic variability this source ex- 
hibits, binning data over multiple nights may flatten the 
lightcurve shap e , impl ying a lower inclination. Similarly, 
iBeekman etafl (|1997l ). who find 10° < i < 26° use a 
binned R-band lightcurve obtained between 1995 Jan- 
uary 12-15. The lightcurve consists of just 55 points, 
binned in phase bins of width 0.1. 

Given the large variability exhibited by GRO 
J0422+32, binning the lightcurves is unusually problem- 
atic. In addition, those measurements fall well below the 
lower limit obtained using the spectral type of the sec- 
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ondary. Therefore we discount measurements obtained 
by binning optical lightcurves. We adopt i = 45°, assume 
the source is active, and adjust the inclination according 
to §4, obtaining an inclination i ~ iV(63.7°, 5.2°). We 
note that the 3tr lower limit o f this distribution is con - 
sistent with the lower limit of iFilippenko et all (|1995h . 
The peak of the probability distribution for the mass of 
the black hole is 2.1M©, which may call into the question 
the nature of this compact object. 

5.5. GRS 1009-45 

The only analysis o f ellipsoidal var i ability in GRS 1009- 
45 was performed bv lShahbaz et all ()1996l ). They obtain 
R-band photometry on 1995 May 8-10 and 1996 Febru- 
ary 17-21. The lightcurve shows evidence of aperiodic 
variabity: the peak brightness fluctuates by around 0.1 
mag during the time of observatio ns, clear evidence of a 
significant nonstellar contribution. IShahbaz et al.l (|1996f ) 
bin the data and fit a star-only model, obtaining a best- 
fit inclination i = 44° with x? e d = 3.2. Such a high y 2 
indicates that the lightcurve is most likely active, as is 
also suggested by the variability. We therefore adjust the 
inclination estimate i = 44° according to §4, obtaining 
i ~ JV(62.0°,5.1°). 

Thi s result conflicts with the work of IFilippenko et all 
(|1999t) . who determine inclination assuming the sec- 
ondary is not undermassive. They find i = 78° for 
a K7-K8 star with mass m* ~ 0.6M Q . According to 
the reasoning in §4.3, this inclination should be a lower 
limit. How ever, there is some amb iguity about the spec- 
tral type. Idella Valle et al.l (|1998l ) suggest the spectral 
type may be as early G5 V. Following the procedure 
outlined in §4, we find that a G5 V star implies a 
lower limit of i = 42.3°, consistent with our estimate 
i ~ JV(62.0°,5.1°). 

5.6. GRS 1124-68 

GRS 1124-68 is one of the few sources for which simul- 
taneou s spectroscopy and photometry exist. lOrosz et al.l 
(|1996| ) analyze spectroscopy obtained 1992 April 3 and 
photometry between 1992 April 3-15 in I and W (the W 
filter is a wide bandpass centered at 4700 A) . Ellipsoidal 
modulation is easily discernible in this data. By con- 
trast, long-term monitoring of the source with SMARTS 
shows no clear ellipsoidal variability since 2003 (C. Bai- 
lyn, private communication). We therefore assume the 
source was passive during the 1992 observations and ac- 
tive post-2003. 

Using the 1992 spectroscopy, lOrosz et al.l (|1996f ) con- 
strain the B+V-band NSL fraction to 0.45 ± 0.05. Be- 
cause the source was passive during this period, the NSL 
fraction is most likely valid for the 12 days of photome- 
try following the spectroscopi c observations. Usi ng NSL 
fractions in the quoted range, lOrosz et al.l (|1996f ) fit the 
1992 B+V lightcurve with a star + constant flux off- 
set model. However, the B+V lightcurve has uneven 
maxima, a signature of the hotspot in the passive state. 
Fits to the larg er maximum yield inclinations above the 
eclipse limit, so lOrosz et al.l (|1996|) fit the smaller maxi- 
mum and obtain 54° < i < 65°. 

This meas urement is somewhat higher and less precise 
than that of IGelino etHI (|2001[ ). They find i = 54 ± 
1.5° from a star-only fit to J and K lightcurves obtained 



2000 February 20-21. However, ther e is clear evidence fo r 
substantial nonstellar flux in the IR IGelino et al.l (|2010D . 
when the source was active. The source may have been 
passive in 2000 February, but that does not guarantee a 
negligi ble NSL frac t ion. 

The lOrosz et al.l (|1996|) result is also slightly higher 
than that of IShahbaz et al.l (|1994l) . who find i = 54^°° 
for a folded, binned H-band lightcurve from 1995 Febru- 
ary 17-20. However, the best fit inclination gave Xred = 
5.4, which the authors suggest is a consequence of in- 
correct sk y subtraction and fiat-fielding. We therefore 
adopt the lOrosz et al.l (|1996f ) result as the most reliable 
inclination measurement, 54° < i < 65°. 

5.7. GRS 1915+105 

GRS 1915+105 is the only X-ray binary for which 
we can mea sure inclination witho ut analysing ellipsoidal 
variability. iGreiner et al.l (120011 ) determine an inclina- 
tion i = 70 + 2° based on the apparent superlumi- 
nal motion of the jet. They assume the jet is per- 
pendicular to the orbital plane due to a lack of ob- 
servable precession over several years of observation. 
GRS 1915+105 has been active since its discovery in 
1994 ([Mirabel fc Rodriguez! 119941 ) . Single- hump vari- 
ability and a su p erhum p period have been detected 
iBailvn fc Buxtonl ([201 If ), but no studies of ellipsoidal 
variability have been possible to date because of sub- 
stantial disk and jet flux variation. 

5.8. GS 1354-64 

ICasares et al.l (|2009l ) obtain multiwavelength photom- 
etry and spectroscopy between 1995 and 2003. The data 
is characterized by strong aperiodic vari ability, and no 
clear e llipsoidal modulation is discernable. ICasares et al.l 
(|2009t) measure a NSL fraction of 0.67 in 2004 and 0.5 in 
2006. The average R-band magnitude varies by nearly 1 
mag be tween 2002 and 2004. Because of this strong vari- 
ability, ICasares et al.l (|2009f) cannot determine a lower 
limit to inclination based on ellipsoidal variations. From 
the spectral type and eclipse limits, we find 27.2° < i < 
80.8° according to the criteria in §4.3. We assume an 
isotropic distribution i ~ 7(27.2°, 80.8°), indicating that 
cosi is uniform between cos 27.2° and cos 80.8°. 

5.9. GS 2000+25 

The most extensive study of e llipsoidal varia b ility i n 
GS 2000+25 was performed by lloannou et all <|2004f ). 
They obtain I-band data from 1992 August and 1998 
September, as well as R-band data from 1998 September- 
October, 1999 June, July, and September, and 2000 July- 
August. The authors note some changes in lightcurve 
shape between the observing runs: for example, the R- 
band lightcurve has equal maxima in 2000 August, but 
unequal maxima in 1998 September. These changes are 
small, however, and there is relatively little scatter in 
the folded lightcurve. We therefore conc l ude th at the 
source is passive during the lloannou et al.1 (|2004[ ) obser- 
vations. The authors bin the lightcurve and fit a star- 
only model, finding 54° < i < 60°. Fits including a disk 
and hotspot give lower Xred values, but do not affect the 
allowed range of inclinations. Fits to data between phase 
0.0 and 0.5, which the authors suggest are less affected 
by the hotspot, also produce the same result. This lack 
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of sensitivity to changes in the model may well be a con- 
sequence of combining and binning data taken over many 
nights: features such as a hotspot are obscured because 
the hotspot position changes o ver time. We suspect that 
binning the data depresses the lloannou et alJ ()2004l ) in- 
clination measurement: even for passive lightcurves, bin- 
ning slightly decreases the amplitude o f the ellipsoidal 
variations. Indeed, the lloannou et al.l (I2004T) measure- 
ment is somewhat lower than that of ICallanan et ail 
(|l996l) . who find 55° < i < 65° for a star-only fit to J and 
K' photometry obt ained on 1995 August 11-12, assum- 
ing q < 0.05. The ICallanan et al.l (|1996l ) result is most 
likely passive, b ecause it gave a highe r inclination range 
tha n the passive lloannou et al.l (120041) lightcu r ves. B oth 
the lloannou et alJ (|2004l ) an d ICallanan et al.j (119961) re- 
sults are consistent with the IBeekman et alJ ( 19961 ) esti- 
mat e 43° < i < 60 ° . For our final estimate, we adopt 
the ICallanan et al.l ()1996l ) result because it is l ess af - 
fected by binning than that of lloannou et al.l (|2004h . 
given that the observations span only two days. We take 
i »j 7(55.0°, 65.0°) as our inclination estimate. 

5.10. GS 2023+338 

Strong aperiodic variability in the lightcurvc of GS 
2023+388 h as thus far prevented precise inclination mea- 
surements. I Wagner et al.l (|1992f ) obtained I-band pho- 
tometry on 17 nights between 1990 September and 1992 
May. During this period, the average I magnitude var- 
ied by more than 0.1 mag, so the authors are only able 
to constrain 50° < i < 80°. The lower limit is derived 
based on the observation of double-peaked Balmer lines 
in the spectrum , and the upp e r limi t is based on the 
lack of eclipses. IShahbaz et al.l (|1994f ) find a similar al- 
lowed range, 45° < i < 83°, using K and K' photome- 
try obtained between 1992 August and 1993 December. 
Their best fit with a star-only model has l = 56.0 with 
a Xred = 15.2. They suggest the poorness of the fit may 
be due to an incorrect color correction between the K 
and K' data. We the refore regard the limits obtained by 
IWagner et ail (|1992l) as more secure 

ISanwal et al.l ( 1996h also fit a star-only model to IR 
data: they obtain H-band photometry on 17 nights be- 
tween 1993 June and November. Strong aperiodic vari- 
ability is present in the lightcurvc. The authors note that 
there is significantly more scatter in the observations of 
GS 2023+388 than in the lightcurvc of the comparison 
star. They also note a systematic increase in brightness 
over one 6-hour period of observation. At the begin- 
ning of that night, the data is fainter than the ellipsoidal 
model, but the source smoothly increases in brightness 
over 6 hours until it is nearly 0.1 mag brighter than the 
model lightcurve. This hour-scale time variability sug- 
gests that the source was active during these observa- 
tions. Such be havior is typical of the a ctive state. On 
the other hand, iKhargharia et all (|2010t ) find almost no 
disk contamination in the IR spectrum. However, these 
data were obtained in 2007, an d thus do not s peak to the 
state of the system when the ISanwal et al.l (|1996h data 
were taken over a decade earlier. We therefore use th e 
minimum inclination estimate of ISanwal et al.l (|1996l ). 
i > 62.0°, which they obtain using a star-only model 
with standard limb and gravity darkening. We assume 
that the source was active at the time. Using the method 
described in §4, we obtain i ~ iV(80.1°, 5.1°). Wc note 



that this is approximately equal to the eclipse limit. 

5.11. H1705-250 

The inclination o f H1705-250 was first measured by 
iMartin et ail (|1995l ). They obtain 87 R-band images on 
1992 May 1-6 and fit a star-only model, finding 48° < 
i < 51°. The best fit had x? ed = °- 87 - The authors 
only show a folded lightcurve, so it is difficult to detect 
a periodic variabili ty. 

iRemillard et al. (19961)) obtain a conflicting result us- 
ing B+V observations from 1992 May, 1993 April, and 
1994 July. The lightcurve exhibits uneven maxima, 
which suggests the source was active during the obser- 
vations. However, the authors suggest this may be an 
artifact of changes in phase coverage between the obser- 
vations, They find a lower limit on i by fitting a star-only 
model to the lightcurve, excluding the higher maximum 
from the fit. Their fit restricts i > 60°. 

Given the uncertainty in whether the source was ac- 
tive or passi v e for both these analyses, we adopt the 
IMartin et ail (|1995l ) i = 48° estimate as a lower limit. 
We cannot obtain an upper limit due to eclipses for 
this source because the limits on the mass ratio extend 
to 0.0. Our final estimate for the inclination is thus 
i ~ 7(48.0°, 90.0°). 

5.12. SAX J '1819.3-2525 



lOrosz et~aTI (|2001l) a nalyze the p hotog raphic B-band 
lightcurve obtained by IGoranskijI (|1990l ). The photo- 
metric errors are large (near 0.1 mag), so the binned 
lightcurve is imprecise. The best-fit model to the 
lightcurve has i = 70° and a partial eclipse of the sec- 
ondary by a large, faint accretion disk; however, it sys- 
tematically underestimates the amplitude of the ellip- 
s oidal variability. 

iMacDonald et al.l (|2011| ) have compiled 10 years of 
data on this source, which they separate into clearly de- 
fined passive and active states. Analysis of the passive 
data shows no nonstellar contribution, a result that is 
supported in some cases by simultaneous spectroscopy. 
However, the amplitude of the ellipsoidal variability re- 
quires an extreme inclination of ~ 90°. The X-ray data 
are sparse enough that an eclipse geometry cannot be 
excluded. Therefore, we adopt i ~ 7(80.0°, 90.0°). 

5.13. XTE J1118+480 

There is consensus in the literature that XTE 
Jl 118+480 has a high inclination, in the range 68° < i < 
81°. However, there are a number of factors that make 
an accurate measurement challenging. One is a strong 
supcrh ump modulation in addition to the ellipsoidal vari- 
ability ([Zurita et al.ll2002l ) . Another comp licating aspect 
is the large and variable NSL fraction: IWagner et alJ 
(|2001[ ) measure an R-band NSL fraction of 0.72 for 2000 
N ovember 30 which d ecreases to 0.64 by 2001 January 4. 

IGehno et ail (|2006l ) obtain data in B, V, R, J, H, and K 
that may be passive. Their R-band lightcurve is a round 
0.8 mag fainter than that of IWagner et al.l (|2001| ). and 
there is no detectable superhump period, both indica- 
tions of less disk activity. They fit data in all bands 
simultaneously and constrain the NSL fraction in each 
band. They find i = 68 + 2°. However, the authors 
assume there is zero disk contamination in the H-band 
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when fitting the spectral template, an assumption that 
may be problem a tic. 

IWagner et "all (|2001h measure i = 81 ± 2° with a 
star+di sk model usin g a brig hter R-band lightcurve than 
that of iGelino et al.l (|2006| ). This is unusual, because 
lightcurves with very high levels of disk contamination 
generally cause lower inclination measurements. How- 
ever, these data were taken on three separate occasions 
between 2000 December 14 and 2001 January 9, dur- 
ing which time the measured R-band disk fraction var- 
ied significantly. Moreover, the lightcurve consists of 68 
exposures that were binned in 25 phase bins. The com- 
bination of binning small amounts of data taken weeks 
apart during a period of intense disk activity may lead 
to an unrepresentat i ve ligh tcurve shape. 

iMcClintock et al.l (12001D also measure a very high in- 
clination, finding a best fit i = 80° for a disk fraction of 
0.66. However, they acknowledge uncertainty in the NSL 
fraction and set a lower limit i > 40° for zero disk light. 

The most extensive study of ellipsoida l varia bility in 
XTE J1118+480 is that of IZurita et "all (|2002l ). They 
obtain R-band photometry on 53 unique nights between 

2000 December and 2001 June. They observe a steady 
decrease in magnitude over this time interval, so they de- 
trend the lightcurve by subtracting the average nightly 
flux. This lightcurve shows evidence of superhump mod- 
ulation and may be distorted by flickering on sub-orbital 
timescales. They infer a NSL fraction for the detrended 
lig htcurve by ext r apolat ing the 2001 April measurement 
of IWagner et all (|2001| ) and assuming the decline in 
brightness is due to a decrease in disk flux. Using the 
extrapolated NSL fraction, they find an inclination in 
the range 71° < i < 82°. 

None of the inclination measurements discussed above 
are free of significant systematic sources of error; they 
are, however, reasonably consistent. We adopt the full 
range of inclination estimates obtained for this object, 
taking i ~ 7(68.0°, 82.0°). 

5.14. XTE J1550-564 

IQrosz et al.l (|2011l) determine the inclination of XTE 
J1550-564 using photometry and spectroscopy obtained 
between 2001 and 2008. The photometric observations 
include optical data from 2001 June and NIR data taken 
between 2006 and 2008. The measured V- and R-band 
NSL fractions are 0.3 and 0.39 using 2001 and 2008 spec- 
troscopy, respectively. There are definite changes in the 
NIR lightcurve shape between 2006-2007 and 2008, so 
it it unlikely that the 2008 NSL fraction measurement 
is valid for the 2006-2007 observations. It is also un- 
likely that the 200 1 NSL fraction is v alid for any of the 
NIR observations. IQrosz et al.l (|2011f ) acknowledge this 
uncertainty, but fit the data using eight different combi- 
nations of lightcurves and NSL fractions. There are four 
data subsets: optical data only, optical data and 2006- 
2007 NIR data, optical data and 2008 NIR data, and all 
optical and NIR data. Each subset is fit separately for 
both NSL fraction measurements. The model includes 
a disk with four free parameters: radius, flaring angle, 
inner radius temperature, and temperature profile. The 
lowest Xrod i s obtained using the optical lightcurves from 

2001 and the 2008 NSL fraction. However, as discussed 
in §3.3, using NSL fractions determined at a different 
time from the observations produces unreliable inclina- 



tion measurements. Nevertheless, the range of inclina- 
tion measurements from all eight combinations is rea- 
sonably narrow: 57.7° < i < 77.1°. Because we do not 
know which of the combinations is most appropriate, we 
adopt i ~ 7(57.7°, 77.1°) as our estimate of inclination. 
We followed the same reasoning to obtain an estimate 
for the mass ratio, assuming q is uniformly distributed 
between the minim um and maximum values found by 
IQrosz et all (|20T1 . 

5.15. XTE Jl 650-500 

IQrosz et al.l (|2004|) determine the inclination of J1650- 
500 with R-band photometry obtained between 2003 
May and August. The authors set a lower limit on 
inclination o f i > 50° by fitting a star-only model. 
IQrosz et al.l (|2004f ) also attempt to constrain the NSL 
fraction using spectroscopy from 2002 June, tentatively 
finding disk contamination near 0.8 in the R-band. There 
is no indication, however, that this measurement is ap- 
plicable to the 2003 photometry. No xf-ed is given for the 
star-only fit, but we infer that the source was active for 
the 2003 observations because there appears to be more 
scatter in the folded lightcurve than one expects due to 
photometric error. We therefore assume the source was 
active and scale the lower limit i = 50° according to §4, 
which implies i ~ N(75.2°, 5.9°). 

5.16. J1859+226 

ICorral-Santana et al.l (|2011| ) obtain R-band photome- 
try on 2008 July 31 - August 1 and 2010 July 13-14. 
T hey compare this d ata to the 2002 R-band lightcurve 
of IZurita et all (|2002l ) and find an increase in brightness 
of ~ 0.25 from 2002 to 2008 and - 1.0 from 2002 to 2010. 
No clear ellipsoidal modulation is detectable in the 2010 
data, but the 2000 and 2008 data show variability with 
amplitude 0.3-0.4 mag, cons istent with the passive state. 
ICorral-Santana et al.l (|2011| ) find that a star-only model 
with i = 60° reproduces these data. We therefore adopt 
this value as representative of the passive state, and es- 
timate i = 60 ±3.0°. 

6. THE BLACK HOLE MASS DISTRIBUTION 

We now analyze the mass distribution of black holes us- 
ing the orbital parameters and new inclination estimates 
listed in Table 2. Previous analyses of the mass distri- 
bution indicate that there is a "mass gap," or dearth of 
black holes, between the maximium theoretical neutron 
star mass (« 3M Q ) and the mimimum black hole mass 
(|Bailvn et al.l \l99§. lOzel et all l2010l: iFarr et alJl20ll . 
The presence of a mass gap has important implications 
for the physics o f black hole formation, as discussed by 
iBelczvnski et al.1 (|2011f ). 

Based on our arguments in §5 that many of these pub- 
lished masses may be overestimates, it is plausible that 
the mass gap inferred in earlier work is the result of sys- 
tematic errors in mass measurements. Using the revised 
system parameters in Table 2, we show in fl6.3l that there 
is no evidence of a mass gap, even when using models of 
the black hole mass distribution that give strong evidence 
for a gap in previous analyses. We demonstrate in i j6.3l 
that this conclusion rests on the properties of one system 
in our sample, GRO J0422+32. However, it is not clear 
that the inclination correction we have identified by anal- 
ogy to A0620-00 is entirely appropriate for this system; 
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TABLE 2 

Orbital parameters for 16 black hole binaries 
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GRS 1915+105 


V1487 Aql 


KO-7 III 


7V(9.5, 3.0) 


7V(0.058, 0.033) 


iV(70.0, 2.0) 


0.0 


90.0 


14,15 


GS 1354-64 


BW Cir 


GO-5 III 


7V(5.73,0.29) 


7V(0.12, 0.04) 


7(27.2,80.8) 


27.2 


80.8 


16,17 


GS 2000+25 


QZ Vul 


K3-6 V 


iV(5.01, 0.12) 


iV(0.042, 0.012) 


7(55.0,65.0) 


28.3 


86.7 


18,19 


GS 2023+338 


V404 Cyg 


K0 IV 


7V(6.08,0.06) 


iV(0.060, 0.005) 


iV(80.1,5.1) 


35.4 


80.0 


20 


H1705-250 


V2107 Oph 


K5±2 V 


7V(4.86,0.13) 


(7(0, 0.053) 


7(48.0,90.0) 


0.0 


90.0 


21,22 


SAX J1819.3-2525 


V4641 Sag 


B9 III 


AT (2. 74, 0.12) 


7V(0.67,0.04) 


7(80.0,90.0) 


44.8 


69.6 


23 


XTE J1118+480 


KV UMa 


K5 V 


iV(6.27,0.04) 


7V(0.024, 0.009) 


7(68.0,82.0) 


21.8 


89.4 


24,25 


XTE J1550-564 


V381 Nor 


K3±l III 


7V(7.65,0.38) 


(7(0.031, 0.037) 


7(57.7, 77.1) 


26.5 


82.0 


26 


XTE J1650-500 




G5-K4 III 


7V(2.73,0.56) 


(7(0,0.5) 


7V(75.2,5.9) 


0.0 


90.0 


27 


XTE J1859+226 


V406 Vul 


K5-7 V 


iV(4.5,0.6) 


(7(0,0.5) 


7V(60.0,3.0) 


0.0 


90.0 


28 



Notes: This table gives the mass function, mass ratio, spectral type, and inclination estimate for 16 black hole binaries. 

The measurements of /, q, and spectral type are taken from the literature. The lower limits on inclination, i m i n , 



highest inclinations that do not cause eclipses. For a detailed explanation of these limits and the inclination 
estimates, see §4. 

Notation: The notation N(fi, a) implies a normal distribution with mean fi and standard deviation cr. A uniform 

distribution from a to /3 is indicated by U(a,/3). An isotropic distribution is denoted i ~ /(o,/3), implying cosi 

is uniform betwe en a and fi. 

References: (DIOrosd pQ03T): (2HOrosz et al.l lH998h; (31 iN cilscn et al] lf2008T): (41 ICantrell et all lIMIoh; ( 51 Shahb aTeTaTI 
(1999); (6) Greene et al. (2001); (7) Webb et al. (2000); (8) Harlaftis et al. (1999); (9) Filippcn ko et alj dl999li : 
flOIIShahbaz et al J 11999); (ll)[d clla Vallc et al. (199J); (1 21IOrosz et all dl996D ; (13) ICasares et all (I1997lr ~(14) 
I Greiner et al.l f2001l); (15) IHarlaftis fc Grein er (2004): (16)ICasares et al.l (1200911 : (17) ICasares et al.l (1200411: ( 18) 
iHarlaftis et alj 119961): (19 1 ICasares et al.l 1199a): ( 20) ICasares fc Charlesl (119941); (21) iFilippenko et al.l 
(221 IHarlaftis et al.l il997|1: (23 ) IOrosz et alj (1200111: (24 TlGonzalcz Hernandez et al. (200811: (25) ICalvelo et"^l 
(120091 1: (261 IOrosz et alj {20111 1: (271 IQrosz et al.l l|2004T l: (28) ICorral-Santana et al.l (|2miT l 



we discuss this issue in M6.3.1I Further observations of 
the system will be needed to settle this question. 

6.1. Statistical Methods 

To address the impact of systematic error in the mass 
distrib ution, we r e peat a subset of the Bayesian anal- 
ysis m iFarr et al.l (|20H using the adjusted inclination 
measurements discussed in §5. Using the distributions 
of system parameters from Table 2, we can compute the 
probability distribution for the true mass of each system. 
Figure 7 shows these distributions for all 16 BHSXTs. 

For a set of parameters, 8, governing the underlying 
BH mass distribution and observational data, D, Bayes' 
rule states 



Pi 



D) 



P(D | 0)P(9) 
PjD) ' 



(8) 



Here P (8 \ D), called the posterior, is the probability of 
obtaining particular values for 8 given the observational 
data, D. P (D \ 8) called the likelihood, is the proba- 
bility of obtaining the observed data, D, when the un- 
derlying mass distribution is described by parameters 8. 
P (8) , called the prior, is the probability distribution one 
would expect for the parameters before seeing the data. 
P(D), called the evidence, is a normalization constant 
that ensures that the posterior is a proper probability 
distribution on parameter space. 

Our data D is the set of probability distributions for 
the underlying masses derived from the observations dis- 
cussed in §5 and displayed in Figure [7] Write Pi (men) 
for the distribution of masses in system i. We assume 



that the mass measurements for the 16 systems are in- 
dependent; then the likelihood is given by 



P(D | 6) = H J F(m BH | e)Pi(m BK )dm B H- (9) 



It is possible to evaluate the posterior in Equation ([8]) 
directly for various values of the parameters, 8. How- 
ever, this calculation is computationally expensive in 
multi-dimensional parameter space. As a more effi- 
cient alternative, we draw parameter samples from the 
posterior distributio n via Markov chain Monte Carlo 
(MCM C) methods (|Metropolis et al.l Il953t IFarr et al.1 
1201 lh . Given a set of samples, {&i}, one can compute 
probability distributions for individual parameters by 
histogram, and can approximate posterior-weighted in- 
tegrals over parameter space by sums, 



(/(f)> P 



WD) 



J d8f(6)p(8\D)^±J2f^- 



(10) 



6.2. Model 

We focus on the power-law model from IFarr et all 
(|2011l Equation (7)) because it was the most-favored 
model for the LMXB mass distribution out of the 10 
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Fig. 7. — Kernel density estimates of the probability distributions for black hole mass implied by the system parameters in Table 2 
(solid lines). The distributions implied by the system parameters used in Fair ct al. (2011) arc also shown (dashed lines). Each panel is 
normalized so that the probability distributions integrate to 1. 



considered in that work 1 . The power law model has 

P («1BH| {A/min, M max , a}) = 

'Amg H M min < m B H < M, 







otherwise 



with parameters 



9 = {A/ min , A/ max ,o;}. 



, (11) 



(12) 



1 We have verified that the qualitative behavior of the mass 
gap described in this section is present for the oth e r mod els of the 
LMXB mass distribution considered in lFarr et al.l 1 1201 II) . 



The normalization constant, A, is 

1 + a 



A 



M 1+a 



M 1+a ' 

mm 



(13) 



We use uniform priors on M m i n , M max > M m - m , and a 
within broad ranges that allow for black hole masses up 
to 40 M Q : 

P(8) = 



9^r— < M 
otherwise 



< M max < 40, -15 < a < 13 



(14) 
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In the analysis of iFarr et all (|2011l ). the power- law 
model had strong evidence of a mass gap between the 
theoretical maximum mass of the heaviest neutron stars 
(~ 3M Q ) and the mass of the lightest black holes, de- 
fined as the 1% q uantile (M-\ %) of t he mass distribution 
in Equation dTTJ. IFarr et alJ (|2011l ) had M 1% > 4.3M 
with 90% confidence for the power-law model. 

6.3. The Mass Gap 

With the revised system parameters in Table 2, there 
are two systems whose mass distribution peaks below 
4M : 4U 1543-47 and GRO J0422+32. Thus it is not 
surprising that the expected mass distribution under the 
power-law model, defined by 

(P) (m BH ) = / d6P(m BH | 6)P{9 \ D) (15) 

gains support in the mass gap. That is, the existence of 
two sources whose mass distribution is strongly peaked 
below AMq necessarily requires that the overall mass dis- 
tribution not be zero in that r egion. In FigureEl we show 
the (P) (mBR-) implied by the lFarr et al.l (|20 111) analysis 
and the (P) (jubh) implied by the system parameters in 
Table 2. Interestingly, the extra support in the gap re- 
gion is due almost completely to the shift in the mass 
distribution of GRO J0422+32; in Figure[8]we also show 
(P) (meii) with parameters from Table 2, but excluding 
GRO J0422+32 from the sample, and similarly for ex- 
cluding 4U 1543-47. When 4U 1543-47 is excluded from 
the analysis, the curve is essentially the same as when 
the entire set of 16 systems is analyzed, while excluding 
GRO . 10422+32 produce s a curve that is very close to 
that of IFarr et all (|2011( ). 

Another way to address the presence of the mass gap 
is to examine the posterior probability distribution of 
the 1% mass qu antile, M-[%. Figu re [HI shows this distri- 
bution from the lFarr et al.l (|2011[ ) parameter values, the 
updated parameters in Table 2, and using the updated 
parameters but excluding 4U 1543-47 or GRO J0422+32. 
The analysis using the complete set of systems from Ta- 
ble 2 has significant probability for Mi% < 3M ; so does 
the analysis with parameters from Table 2, but exclud- 
ing 4U 1543-47. On the other ha nd, the analysis with 
parameters from lFarr et al.1 ([201 ID and Table 2, but ex- 
cluding GRO J0422+32, have M x % > 4.3M with 90% 
confidence. Thus, the presence or absence of a mass gap 
depends strongly on the properties of the mass distribu- 
tion for GRO J0422+32. In addition, if in the future the 
mass distribution of 4U 1543-47 were constrained to a 
narrow peak near the low-mass end of its current distri- 
bution, this system would also provide strong evidence 
against a mass gap. Because the continued evidence for 
or against the mass gap may depend on the constraints 
on these systems from future observations, we next dis- 
cuss the current status of observations of these systems 
in more detail. 

6.3.1. GRO J0422+S2 

The inclination of GRO J0422+32 hinges on whether 
this source is generally observed in an active state, and 
whether that active state biases the inclination measure- 
ments in a similar way to that of A0620-00. As shown 
above, if the typical measurements of i ~ 45° are ac- 
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Fig. 8. — The marginalized mass distribution, Equation 1151 1. for 
the power-law mass distribution model (Equation {TTJ). The solid 
curve from an analysis using the system parameters in Table 2; the 
dotted curve from an analysis with the parameters from Table 2, 
but excluding GRO J0422+32 from the sample; the dashed curve 
comes from an analysis with system parameters from Table 2, but 
excluding 4U 1543-47 from the sample; and the dot-da shed curve 
comes from the analysis using the system parameters of lFarr et al.l 
II2011D ; The analysis using parameters from Table 2 shows sig- 
nificant support in the gap region, but this conclusion depends 
strongly on the updated parameters for GRO J0422+32. 

cepted, then the previously identified mass gap still ex- 
ists. However, if it is larger, as is suggested by the anal- 
ogy with A0620-00, then the compact object in this sys- 
tem must be roughly between 2 and 3 solar masses, which 
effectively fills the mass gap (see the relevant panel in 
Figure 7). In Figure [T0l we show how the 16% quantile 
of J0422+32's mass distribution varies as its inclination 
distribution's mean varies from the previously-accepted 
value to the value in Table 2, holding the width of the 
inclination distribution fixed. The 16% quantile approx- 
imately tracks the peak of the M\% distribution shown 
in Figure [3J so we can see the effect on the mass gap 
from varying the inclination without repeating the anal- 
ysis for each value of i. We expect to obtain no evidence 
of a mass gap as long as the inclination distribution of 
J0422+32 peaks at i > 45°. 

In many ways, GRO J0422+32 is just the kind of sys- 
tem one would expect to be heavily biased by non-stellar 
light. The secondary star is one of the smallest and 
coolest of the entire sample. Therefore changes in the 
non-stellar flux will be reflected strongly in the over- 
all flux of the system. And indeed, GRO J0422+32 is 
quite faint in quiescence (i? ~ 21) and strongly vary- 
ing lightcurves are observed at different tim es (com- 
pare [Ounxii^ Thus it 
seems likely that the observed lightcurves of this source 
are distorted by non-stellar light in ways that change 
strongly with time. Whether these distortions conspire 
to create significant underestimates in the inclination, as 
in A0620-00, is not yet clear. One route to resolving 
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Fig. 9. — The distribution of the 1% mass quantile, M]%, in the 
power-law model implied by the analysis with system parameters 
from Table 2 (solid line); parameters from Table 2, but excluding 
GRO J0422+32 (dotted line); parameters from Table 2, but ex- 
cludin g 4U 1543-47 (dashed line); and parameters from lFarr et alj 
1)20111 ) (dot-dashed line). For the complete set of parameters from 
Table 2, and this set excluding 4U 1543-47, the range of likely min- 
imum black hole masses extends through the gap (M195 < 4Mq); 
when GRO J0422+32 is exclud ed the min i mum black hole mass is 
equivalent to the analysis from Farr et al. (2011). 
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Fig. 10.— The variation in the 16% quantile of J0422+32's 
mass distribution, i7i 16 % , with the assumed peak of the inclination 
distribution for this object. m 16 % is a rough proxy for the peak of 
the Mi% distribution (see Figure Q5). As long as the inclination of 
J0422+32 is peaks at i > 45° we expect to obtain no evidence for 
a mass gap. 



this issue would be to carry out simultaneous photome- 
try and spectroscopy of GRO J0422+32, so that the di- 
vision of stellar and non-stellar light can be determined 
for each individual observation, resulting in an unbiased 
light curve of stellar light only. However, given the faint- 
ness of this source, such observations will be challenging. 

6.3.2. 4U 1543-47 

The other source which has a mass probability distribu- 
tion significantly lower than in previous wor k is 4U1543- 
47. T his is because we did not accept the lOrosz et al.l 
(|2002D inclination value, since no published lightcurve 
accompanied that report, but rather accepted the ear- 
lier, significantly higher estimate of Orosz et al. (1998). 
That previous higher inclination measurement resulted 
in a significantly smaller black hole mass, but still al- 
lowed masses above 4M Q , so the source does not greatly 
affect the possibility of a mass gap. However, if the ob- 
servations become more precise at the higher inclination 
value this source might also be an issue for the mass gap, 
so we discuss it further here. 

The difficulties with 4U1543-47 are associated primar- 
ily with its low inclination. The sin 3 i term in equation 
(1) means that at low inclination very small changes in i 
cre ate large cha n ges in the black hole mass. Going from 
the lOrosz et all (|2002f ) result of i ~ 20° to the previous 
result centered around i ~ 30° results in a change in the 
centroid of the probability distribution of the black hole 
mass of more than a factor of two. At low inclinations, 
the amplitude of the ellipsoidal modulation is also quite 
small (< 0.1 magnitudes in this case), so high precision 
magnitude measurements are required. Given that er- 
rors of a few degrees in inclination determinations seem 
inevitable, it may prove difficult to pin down the inclina- 
tion of 4U 1543-47 sufficiently tightly to fully resolve the 
issue. 

4U 1543-47 does have some observational advantages, 
however. It is one of the three sources (with SAX 1819.3- 
2525 and GRO J1655-40) that have early-type secon- 
daries. These systems all have very regular ellipsoidal 
variations in the passive state, with very little non-stellar 
light; this is presumably because the relatively hot sec- 
ondary stars are particularly bright. Like the other early 
secondary systems, 4U1543-47 is quite bright in quies- 
cence (V ~ 15) so that also makes the system easy 
to study. There is one technical difficulty in studying 
4U1543-47: it is in the wings of an even brighter star, 
but this can be dealt with through appropriate use of 
PSF fitting techniques. A long term photometric data 
set for this source, si milar to those for A0620 -00 (C10) 
and SAX1819.3-2525 (|MacDonald et all 120111) has been 
obtained, and is currently being analyzed. 

6.4. Discussion 

A resolution of the specific situation of GRO J0422+32 
and 4U1543-47 will be required to determine whether a 
true mass gap exists in the BHSXT mass distribution. 
However, even if these two objects do prove to have 
Al < 4Af Q , it is still true that the mass distribution 
rises from the low mass end. Thus the relative paucity 
of low mass black holes in these systems still requires an 
explanation. The peak of the mass distribution appears 
to be around 8M Q , slightly higher than in the previous 
analysis by Farr et al. (2011). The new analysis also 
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confirms another aspect of the previous analysis, namely 
a sharp cutoff above IOMq. This is quite different for the 
BHSXT black hole distribution than for the wind-fed sys- 
tems with M star > M B h (eg. Cyg X-l and M33 X-7), 
which often contain much more massive black holes. The 
difference presumably arises due to the dramatically dif- 
ferent evol utionary scenarios that create these two kinds 
of systems (|Valsecchi et al.ll2010HWong et al.ll2012t ). It is 
important to avoid combining these two kinds of systems 
into a single mass distribution, since there are clearly sig- 
nificant empirical and evolutionary differences between 
them. 

The characteristics of the BHSXT mass distribution 
should provide guidance and constraints for the physics 
of massive core collapse. Previous studies have dis- 
cussed possibilities for observational biases and/or physi- 
cal mechanisms that can explain the charact eristics of the 
low-mass BH distribution in these systems dBrown et al.l 
l200lHFryer_fc Kalogerall200lHOzel et"alll2010t ). Most re- 
cently, iBelczvnski et al.l (|2011h have carefully assessed 



how this low-mass tail and the potential presence of a 
gap leads to constraints regarding the timescale for the 
instability growth that eventually drives the stellar ex- 
plosion, as well as the role of binary evolution in affected 
the BH masses in low-mass X-ray transients. Our current 
results, with the mass gap partially filled in, still provide 
an empirical basis for exploring the physical processes 
considered in these studies, although the constraints on 
the explosion timescale will be weaker. 

7. CONCLUSIONS 

We have explored the systematic effects of non-stellar 
light on the mass estimates of black holes in BHSXTs. 
We summarize our conclusions as follows: 



inclinations, and thus accurate black hole masses, 
we identify a number of problematic issues. 

These problems stem from the presence of variable 
non-stellar light, even in the IR. The amount of 
this extra flux can vary significantly from one ob- 
servation run to another, and thus it is not safe to 
use a measurement of stellar light fraction from one 
time to calibrate photometric observations made at 
another time. The non-stellar light also appears to 
vary with orbital phase, so the assumption that it 
is a flat "dilution" of the ellipsoidal variations is 
not always true. Significant variability also means 
that phase binning over long periods of time may 
result in a lightcurve that is unrepresentative of the 
underlying ellipsoidal variations. 



4. When we apply a consistent set of criteria to the 
existing literature on black hole SXTs, we find that 
two objects, GRO J0422+32 and 4U1543-47, may 
have small black hole masses (below 4-5 Mq, see 
Figure 7). Specifically, the low mass required for 
GRO J0422+32 eliminates the mass gap identified 
in previous work on the black hole distribution. 
But this result depends on the active state of GRO 
J0422+32 biasing the inclination measurement in 
the same way as in A0620-00. While this seems 
plausible, additional observational work will be re- 
quired to assess whether the active state is indeed 
biasing the black hole mass in this system. We note 
that if GRO J0422+32 is excluded from the analy- 
sis presented here, the results from previous studies 
are reproduced quite closely. 



1. By examining the case of A0620-00, for which the 
most extensive and carefully analyzed data set ex- 
ists, we find that observations in the "active" state 
significantly underestimate the inclination of the 
system, and consequently overestimate the black 
hole mass. By contrast, observations in the "pas- 
sive" state appear to be relatively unbiased toward 
the system inclination. 

2. We estimate how large such effects might be in 
other systems, based on the assumption that the 
accretion flows are roughly similar in form, while 
the temperature of the star and the relative size of 
the star and the accretion disk vary from system to 
system. The assumption of similar accretion flows 
appears to be in good agreement with spectroscopic 
measurements of the NSL fractions. 

3. We re-examine the literature on the sixteen known 
black holes in SXTs. The data and analysis of 
many of these systems has been aimed at establish- 
ing a firm minimum for the mass of the compact 
object, since that is what establishes the identity 
of the compact object RS R black hole. When exam- 
ined from the point of view of determining accurate 



5. However, even if GRO J0422+32 (and possibly 
4U1543-47) do fall into the purported mass gap, 
the basic features of the mass distribution of BH- 
SXTs remain: there are relatively few low mass 
(< 5Mq) black holes; there is a peak in the dis- 
tribution around 7 — 8Af Q ; and a sharp dropoff in 
numbers beyond 10M Q . This distribution is quite 
different from the "high-mass" black hole binary 
systems, and provides interesting constraints on 
the supcrnovae and binary evolution processes that 
create BHSXTs. 
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